Infrared phonon anomaly and magnetic excitations in single-crystal Cu 3 Bi(Se03)202Cl 
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Infrared reflection and transmission as a function of temperature have been measured on single 
crystals of Cu3Bi(Se03)2C>2Cl. The complex dielectric function and optical properties along all 
three principal axes of the orthorhombic cell were obtained via Kramers-Kronig analysis and by 
fits to a Drude-Lorentz model. Below 115 K, 16 additional modes (8(E|| a)+6(E|| 6)+2(E|| cj) 
appear in the phonon spectra; however, powder x-ray diffraction measurements do not detect a new 
structure at 85 K. Potential explanations for the new phonon modes are discussed. Transmission in 
the far infrared as a function of temperature has revealed magnetic excitations originating below the 
magnetic ordering temperature (T c ~24 K). The origin of the excitations in the magnetically ordered 
state will be discussed in terms of their response to different polarizations of incident light, behavior 
in externally-applied magnetic fields, and the anisotropic magnetic properties of Cu3Bi(SeOa)202Cl 
as determined by d.c. susceptibility measurements. 



I. INTRODUCTION 

Geometrically frustrated materials possessing mag- 
netic order have garnered interest from both scientific 
and technological viewpoints due to their multiplicity 
of low energy nearly-equivalent states, which in turn 
prevent simple long-range magnetic order (e.g., antifer- 
romagnetism). More specifically, because neighboring 
antiferromagnetically-ordered moments in a frustrated 
structure cannot simultaneously minimize their local ex- 
change energies, they are forced to cant with respect 
to their preferred directions. The canting of magnetic 
moments in a frustrated system is known to produce 
noncollinear short-range magnetic orders^ (e.g., cycloidal 
or spiral spin orders) that in turn yield many device- 
applicable phenomena such as magnetically driven ferro- 
electricity and large magnetoelectric responses ^ In ad- 
dition, the ability to synthesize geometrically frustrated 
magnetic materials in a layered structure facilitates the 
study of the intriguing phenomena associated with frus- 
tration because it allows for a quasi two-dimensional 
model of a three-dimensional system. 

Cu3Bi(Se03)202Cl is a geometrically-frustrated lay- 
ered material possessing magnetic order. P. Millet et alA 
determined the room temperature crystal structure of 
Cu 3 Bi(SeC>3)202Cl using single crystal x-ray diffraction. 
The Cu3Bi(SeC>3)202Cl compound crystallizes in a lay- 
ered structure where the individual layers stack along the 
c axis of the orthorhombic cell (Pmmn). The magnetic 
Cu 2+ ions form a hexagonal arrangement within each 
plane that is reminiscent of a Kagome lattice, and thus 
implies magnetic frustration. Within any one hexagon 
there exist two unique copper sites, Cul and Cu2, which 
possess different out-of-plane oxygen bonding. Both cop- 
per sites form distinct [CUO4] units that are linked by 
Se 4+ and Bi 3+ ions. The Cl~ atoms, which sandwich the 
planes formed by the copper hexagons, position them- 



selves along the axis that defines the parallel stacking 
of the copper hexagons. Pictures and more detailed de- 
scriptions of the crystal structure are found elsewhere^ 
P. Millet et alA also reported magnetic sus- 
ceptibility measurements on a powder sample of 
Cu3Bi(Se03)202Cl. Near 150 K they observed a change 
in slope of 1/%(T) that resulted in two distinct and pos- 
itive Weiss temperatures. To reconcile this anomaly in 
l/x(T), they performed linear birefringence on a single 
crystal and deduced that a second-order structural tran- 
sition was likely, but they were not able to determine 
the exact nature of the transition. P. Millet et al. also 
reported ferromagnetic-like behavior below T c « 24 K. 

The technique of infrared spectroscopy lends itself well 
to the investigation of geometrically frustrated magnetic 
materials. For example, the spin-driven Jahn- Teller ef- 
fect in the Cr spinel compounds CdCr204 and ZnCr204 
acts to lift the degenerate ground state arising from com- 
peting magnetic interactions via a coupling to the lat- 
tice degrees of freedom; the ensuing lattice distortion 
is unambiguously observed in the infrared as a split- 
ting of infrared-active phonon modesi&£ Furthermore, 
the a.c. electric (magnetic) field of the infrared light can 
interact with ordered moments and excite an electro- 
magnon (magnon) excitation, which in turn can provide 
information about the symmetry and nature of the mag- 
netic order £2 

Here we present our infrared studies on single crystal 
Cu 3 Bi(Se03)202Cl. We observe 16 new modes in the 
phonon spectra originating below 115 K. Strikingly, our 
subsequent powder x-ray diffraction measurements reveal 
the same 300 K structure existing at 85 K. Preliminary 
Raman measurements suggest that a loss of inversion 
symmetry is likely. In addition, we observe new excita- 
tions in the infrared arising in the magnetically ordered 
state (below 24 K) that show isotropic infrared polariza- 
tion dependence but anisotropic external magnetic field 
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dependence. In light of the novel excitations observed in 
the magnetically ordered state, we performed d.c. suscep- 
tibility measurements to examine the anisotropic mag- 
netic properties of Cu3Bi(Se03)202Cl. The results of 
our magnetic susceptibility measurements are in agree- 
ment with a recent report*^ of the magnetic properties in 
the similar Cu3Bi(Se0 3 )20 2 Br compound (T c = 27.4 K). 
(This manuscript is accompanied by supplementary ma- 
terial including additional figures referred to in the text 
and the results of full Rietveld refinements.) 



II. EXPERIMENTAL PROCEDURES 

Single crystals of Cu 3 Bi(Se0 3 ) 2 2 Cl were grown by 
standard chemical vapor-phase method. Mixtures of ana- 
lytical grade purity CuO, Se0 2 and BiOCl powder in mo- 
lar ratio 3:2:1 were sealed in quartz tubes with electronic 
grade HC1 as the transport gas for the crystal growth. 
The ampoules were then placed horizontally into a tubu- 
lar two-zone furnace and heated at 50°C/h to 450° C. 
The optimum temperatures at the source and deposition 
zones for the growth of single crystals were found to be 
480°C and 400°C, respectively. After four weeks, many 
tabular green Cu3Bi(Se03) 2 2 Cl crystals with a maxi- 
mum size of 15 x 12 x 1 mm 3 were obtained, which were 
indcntified as synthetic francisite on the basis of x-ray 
powder diffraction data. 

Zero-field temperature-dependent (7-300 K) re- 
flectance and transmittance measurements were collected 
on a 6.7 x 4.8 x 0.7 mm 3 single crystal (crystal 1) using 
a Bruker 113v Fourier transform interferometer in con- 
junction with a 4.2 K silicon bolometer detector in the 
spectral range 25-700 cm -1 and a nitrogen cooled MCT 
detector from 700-5,000 cm -1 . Room temperature mea- 
surements from 5,000-33,000 cm -1 were obtained with a 
Zeiss microscope photometer. Magnetic field-dependent 
transmission measurements in the spectral range 15-100 
cm -1 were performed on a 7.5 x 3.8 x 0.2 mm 3 single crys- 
tal (crystal 2) at beam line U4IR of the National Syn- 
chrotron Light Source, Brookhaven National Laboratory, 
utilizing a Bruker IFS 66-v/S spectrometer. The crystal 
was placed in a 10 T Oxford superconducting magnet and 
the transmitted intensities were measured using a 1.8 K 
silicon bolometer detector. Additional magnetic field- 
dependent transmission measurements were performed 
on a 5.1 x 3.5 x 0.3 mm 3 single crystal (crystal 3) in the 
spectral range 10-45 cm -1 using a modified Polytec FIR 
25 spectrometer interfaced to a 7 T Oxford split-coil su- 
perconducting magnet at the University of Wollongong. 
The transmitted intensities were subsequently converted 
to absorption coefficients using a modified Beer-Lambert 
law that accounted for reflection losses at the surfaces. 
Linearly polarized light was used in all infrared measure- 
ments and oriented along the three principal dielectric 
axes of the system. For an orthorhombic system, the 
three principal dielectric axes coincide with the three 
crystallographic axesJ^ In the following report, the sym- 



bol "E" always refers to the polarization of the incoming 
light, while "H" is the orientation of the external field 
(when applicable). 

Powder x-ray diffraction measurements were per- 
formed at beamline X16C of the National Synchrotron 
Light Source, Brookhaven National Laboratory. A small 
crystal was crushed in a mortar and pestle, mixed with a 
small amount of Si powder (NIST Standard Reference 
Material 640c) and roughly three times its volume of 
ground cork, and loaded in a thin-walled glass capillary 
of 1 mm nominal diameter. The cork served to dilute 
the sample so that it could fill the cross-section of the 
x-ray beam without drastically absorbing it; measured 
transmission at the center of the capillary was 12%. X- 
rays of nominal wavelength 0.6057 A were selected by 
a channel-cut Si(lll) monochromator before the sam- 
ple; small drifts in the x-ray wavelength were corrected 
via the Si internal standard. The diffracted beam was 
analyzed by a Ge(lll) crystal and detected by a com- 
mercial Nal(Tl) detector. Data were typically collected 
in steps of diffraction angle 28 of 0.005°. For lattice pa- 
rameter measurements as a function of temperature, the 
sample was inside a Be heat shield in a closed cycle He 
refrigerator, which was rocked several degrees at each 28 
step. Data for crystallographic refinements were collected 
with the sample continuously spinning several revolutions 
for each point, in an Oxford Cryostream sample cooler. 
Powder x-ray diffraction data were analyzed using Topas- 
Academic software*^ 

Magnetic measurements were performed in a commer- 
cial superconducting quantum interference device mag- 
netometer on crystal 2. 



III. RESULTS AND ANALYSIS 

A. Zero Field Reflectance and Transmittance 
Spectra 

The temperature-dependent reflectance spectra of 
Cu 3 Bi(Se03) 2 2 Cl along the a, b, and c axes in the fre- 
quency interval 30-1,000 cm -1 (4-120 meV) is shown in 
Fig. [TJ Remarkably, many new phonon modes are ob- 
served in the reflectance spectra along all three crystal 
directions upon cooling from 120 K to 110 K. The ar- 
rows in Fig. [T] indicate the position where new modes 
occur. Meticulous temperature sweeps (every 1 K) be- 
tween 120 K and 110 K indicate that all the new modes 
arise at 115 K; this trend was verified upon cooling and 
warming. 

In addition, a strong sharpening of many modes is ob- 
served with decreasing temperature. In regions where 
new modes appear, the typical softening of resonance 
frequencies with increasing temperature is not system- 
atically observed. The typical softening of resonance fre- 
quencies with increasing temperature, which is correlated 
to the expansion of the lattice, can be counteracted by 
repulsion due to phonon mixing. An explanation of the 
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repulsion of phonons, which are bosonic excitations, will 
be presented in the discussion section. In addition to 
the new phonons, we observe the softening (at 5 K to 
^30% of the 300 K value) of the lowest frequency mode 
along the b direction, an occurrence which is not asso- 
ciated with phonon repulsion. The soft mode behavior 
exhibited by this low frequency phonon is reminiscent of 
that occurring in a displacive ferroelectric material. 

Cu3Bi(Se03)202Cl is opaque between ~40 and 
800 cm -1 due to the strong optical phonon absorptions in 
this region; however, above the optical absorptions light 
is observed to transmit. Figure [5] displays the zero field 
300 K transmission spectra of Cu3Bi(Se03)202Cl along 
the a and b axes in the mid-infrared and near-infrared 
regions. (Due to the geometry of the sample no light was 
observed to transmit along the c axis.) The anisotropy, 
which is prominent in the infrared reflectance spectra, is 
also observed in transmission at high frequencies. Sharp 
but weak absorption features exist at 1950 cm -1 along 
the a axis and at 2000 and 2050 cm -1 along the b 
axis. The aforementioned absorptions exhibit a minimal 
strengthening with decreasing temperature. The down- 
turn in transmission at 9,000 cm -1 (~1.1 eV) indicates 
the onset of electronic absorptions. The energy of the 
observed gap is consistent with the insulating nature of 
the material. Transmission below the optical absorptions 
(below 40 cm -1 ) is the subject of section III E. 



B. Kramers-Kronig and Oscillator- model fits 

Along the a and b axes we have utilized a combined 
reflection and transmission analysis to extract the sin- 
gle bounce reflectance necessary for the Kramers-Kronig 
transformation. More specifically, in the frequency re- 
gion where the crystal transmits along the a and b axes, 
the measured reflectance has an additional contribution 
from the back surface. Using the method employed by 
Zibold et al.^ i.e., assuming that k is small in the region 
of interest, we are able to extract the single bounce re- 
flectance of the sample. Since the c axis does not exhibit 
transmission, the measured reflectance was assumed to 
be single bounce reflectance along this direction. 

The real and imaginary parts of the dielectric func- 
tion were estimated from the single bounce reflectance, 
R s (co), using the Kramers-Kronig transformation^ Be- 
fore calculating the Kramers-Kronig integral, the low fre- 
quency (0.1-30 cm -1 ) data were approximated using the 
dielectric function determined from the fitting procedure 
described below. At high frequencies the reflectance was 
assumed to be constant up to 1 x 10 7 cm -1 , after which 
R ~ (w)~ 4 was assumed as the appropriate behavior 
for free carriers. The optical properties were obtained 
from the measured reflectance and the Kramers-Kronig 
derived phase shift on reflection. 

The single bounce reflectance was fit with a Lorentz 
oscillator model to obtain a second estimate of the com- 
plex dielectric function in the infrared range. The model 



assigns a Lorentzian oscillator to each phonon mode in 
the spectrum plus a high frequency permittivity, f^, to 
address the contribution of electronic absorptions. The 
model has the following mathematical form: 

00 S uj 2 

ojj 2 - uj 2 - lujjj 

where Sj , u)j , and jj signify the oscillator strength, cen- 
ter frequency, and the full width at half max (FWHM) of 
the j th Lorentzian oscillator. The reflectivity is then cal- 
culated using the Drude-Lorentz complex dielectric func- 
tion. A figure comparing the calculated and measured 
reflectivity is available in the supplementary information. 

C. Optical Properties 

The optical properties extracted from both Kramers- 
Kronig analysis and our fits of reflectance can be used to 
investigate further the appearance of new phonon modes 
in the infrared spectra at 115 K. The upper panel of 
Fig. [3] depicts the real part of the optical conductivity, 
cti(w), along the b axis in the frequency range 250-285 
cm -1 . The lower panel of Fig. [3] depicts the loss function, 
Im(-l/e(w)), in the same region. Both optical properties 
shown are from Kramers-Kronig analysis of the single 
bounce reflectance. The peaks observed in u\{uj) and 
Im(-l/e(u;)) closely correspond respectively to the TO 
and LO phonon frequencies. The inset of the lower panel 
of Fig.[3]depicts the plasma frequency (fl s = \J SjUj 2 ) as- 
sociated with both phonons as a function of temperature 
as obtained by our Drude-Lorentz fitting. The equation 
defining the plasma frequency of ionic vibrations is sim- 
ilar in form to that which describes the free carrier re- 
sponse in metallic systems after the electronic mass and 
charge are replaced with the reduced mass of the normal 
mode and ionic charge. The plasma frequency is of par- 
ticular interest here because its square is proportional to 
the spectral weight associated with a particular phonon. 
The new phonon mode that appears at ~ 276 cm -1 below 
115 K gains spectral weight at the expense of the existing 
mode at ~ 256 cm -1 . The smooth shift of spectral weight 
with decreasing temperature from the existing mode to 
the new mode is reminiscent of a second order transition. 



D. Powder X-Ray Diffraction 

The changes in the infrared spectra, namely the ap- 
pearance of new vibrational modes below 115 K, sug- 
gest a reduction in lattice symmetry. To investigate the 
low temperature structure of Cu3Bi(Se03)202Cl, we per- 
formed powder x-ray diffraction measurements between 
30 and 300 K. The lattice parameters as a function of 
temperature extracted from Rietveld fits of the diffrac- 
tion spectra are shown in Fig. SI There exists some struc- 
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FIG. 1. (Color online) The 
temperature-dependent 
reflectance spectra of 
Cu 3 Bi(Se03)20 2 Cl along 
the a, b, and c axes. Crimson 
arrows indicated the positions 
1000 Q f new phonons arising at 
115 K. 



Energy (eV) 




1000 



Frequency (cm 1 ) 



10000 



FIG. 2. (Color online) Transmission along the a and b axes of 
crystal 1 in the mid-infrared and near-infrared regions. The 
inset is a sketch of crystal 1 with the a, b, and c axes indicated. 



ture in the curves of lattice parameters vs. temperature, 
but there is nothing that could be regarded as conclusive 
evidence of a phase transition. The negative thermal ex- 
pansion observed for the a lattice parameter below 100 K 
is not unusual for an ionic compound of orthorhombic 
symmetry where it may be energetically favorable for the 
unit cell to contract in one direction while expanding in 
other directions. 

A figure depicting Rietveld refinements of the powder 
x-ray diffraction pattern at 295 and 85 K is available in 
the supplementary information. Both refinements were 
consistent with the published structure of P. Millet et 
al. in space group Pmmn. Full details of the refine- 
ments, including bonding geometry, may also be found 
in the supplementary information. We saw no evidence 
for splitting or broadening of diffraction lines nor the ap- 
pearance of new diffraction lines at low temperature. We 
conclude that there is no direct evidence for a lowering of 
crystallographic symmetry near 115 K. This is a serious 
puzzle because the appearance of 16 new infrared modes 
below 115 K implies that the symmetry of the crystal is 
lower than Pmmn below that temperature. 
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FIG. 3. (Color online) The temperature dependent opti- 
cal conductivity ai(uu) (upper panel) and loss function Im(- 
l/e(io)) (lower panel) along the b axis in the frequency range 
250-285 cm -1 . The insets of the upper and lower panels de- 
pict respectively the temperature dependence of the resonance 
and plasma frequencies of the two modes. 



E. Magnetic Field-Dependent Transmission 

Far-infrared transmission as a function of temperature 
and external magnetic field was measured on crystal 2 
and crystal 3 with light polarized along the a axis, b axis, 
and at 45° to both the a and b axes. Upon cooling to 
5 K an excitation was observed at 33.1 cm -1 that only 
existed in the magnetically ordered state (below 24 K). 
The excitation was observed in all four polarizations. A 
figure depicting the excitation, at 5 K, in each of the four 
polarizations measured, as well as the evolution of the 
excitation as the external field is ramped to 10 T (H \\ c 
geometry) is available in the supplementary information. 

When the external magnetic field was applied parallel 
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FIG. 4. (Color online) Lattice parameters and unit cell vol- 
ume vs. temperature. Error bars reflect the reproducibility of 
independently measured data sets at selected temperatures. 
Solid lines are drawn as guides to the eye. Lattice parameters 
at 30 K are a — 6.3463(2) A, b = 9.6277(3) A, c = 7.2186(3) 
A, Volume = 441.10(2) A 3 . 



to the c axis (H || c) and the field was ramped to 10 T, 
the excitation at 33.1 cm -1 increased linearly with in- 
creasing field (inset Fig. [5^) • No hysteresis effects were 
observed upon ramping the magnetic field. To establish 
confidence that the observed excitation only existed in 
the magnetically ordered state, the crystal was warmed 
up to 40 K and the external field was ramped from to 
10 T. The excitation was not observed. The results are 
shown in Fig. [5£i for the E|| b + 45° polarization. 

In addition, for fields 1 T and greater applied paral- 
lel to the c axis, a second magnetic excitation was ob- 
served which also possessed isotropic polarization depen- 
dence in the ab plane. The excitation increased linearly 
from 10.5 cm" 1 at 1 T to 19.3 cm" 1 at 10 T. Because 
of the low signal to noise level it is not known whether 
the excitation disappeared below 1 T or it was just un- 
resolved. The excitation is pictured for the E \\ a polar- 
ization in Fig. |51 The spectra in Fig. [S] is a conjunction 
of data from two different interferometers on two sep- 
arate crystals of Cu3Bi(Se03)202Cl (see Experimental 
Procedures section). The complimentary techniques on 
separate crystals establish confidence in the existence of 
the excitation. 

When the external magnetic field was applied perpen- 
dicular to the c axis (H _L c) and ramped to 10 T, the 
excitation at 33.1 cm" 1 decreased its resonance frequency 
quadratically with field (inset Fig.[5jD). The field depen- 
dence of the excitation for the E|| a + 45° polarization 
is depicted in Fig. [3b. It should be noted that the de- 
tails of our experimental setup required that the external 
magnetic field be applied parallel to the polarization of 
the incoming light in the H _L c geometry. No additional 
modes in this orientation of external field were detected. 



FIG. 5. (Color online) Experimental support (a) for claiming 
that the magnetic resonance observed at 33.1 cm -1 exists be- 
low the magnetic ordering temperature (T c ~24 K). External 
fields, which act to sharpen the resonance, were also applied 
above T c for further verification that the excitation was not 
observed. In the inset, the dashed red line denoted Hmm indi- 
cates the region where the metamagnetic transition occurs in 
the H|| c geometry. The isotropic magnetic excitation's field 
dependence in the // 1 c geometry (b) with light polarized 
along the E|| a + 45° direction. 



Similar experiments in magnetic fields were carried out 
in reflectance geometry; however, the excitations were 
much too weak to give rise to reflectance bands. 
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FIG. 6. (Color online) The lower frequency magnetic reso- 
nance mode observed in the H|| c geometry for fields of 1 T 
and greater. The two panels result from complimentary tech- 
niques on two separate crystals of Cu3Bi(Se03)202Cl. Ex- 
perimental limitations prevented a study of the mode below 
1 T fields. 



F. Magnetic Properties 

The anisotropic response to external magnetic fields 
of the magnetic excitation observed at 33.1 cm -1 
in transmission has inspired an investigation of the 
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anisotropic magnetic properties of Cu3Bi(Se03)2C>2Cl 
through d.c. susceptibility measurements. The results 
of our d.c. susceptibility measurements are depicted in 
Fig. [7] Isothermal magnetization measurements taken at 
5 K exhibit the strong anisotropic response to the direc- 
tion of an external field previously noted in transmission 
measurements. As shown in Fig. [7^,, with H _L c, the re- 
sulting magnetization vs. field loop resembles that of an 
antiferromagnet. On the contrary, with H || c (Fig. [JJd), 
antifcrromagnetic behavior is observed from to 0.1 T, 
followed by a metamagnetic transition from 0.1 to 0.8 T, 
and then ferro or ferrimagnetic behavior persists from 
0.8 to 5 T. Figure [Jj; is an enlarged view of the metam- 
agnetic transition in the H \\ c geometry. The hysteresis 
observed upon sweeping the field is likely associated with 
the metamagnetic transition. To verify further the anti- 
ferromagnetic to ferro or ferrimagnetic transition occur- 
ring in the H \\ c geometry, magnetization was measured 
as a function of temperature in fields of 0.01 and 1 T. The 
results are shown in Fig. [7Ji. The low temperature can- 
cellation of oppositely aligned moments expected for an 
antiferromagnet is observed at 0.01 T, whereas a strong 
ferro or ferrimagnetic magnetization is observed at low 
temperature for the 1 T magnetization data. 
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FIG. 7. (Color online) Isothermal magnetization measure- 
ments at 5 K, which is well within the magnetically ordered 
state, for H _L c (a) and H \\ c (b,c). The magnetization as a 
function of temperature (d) for H || c at fields above and be- 
low the metamagnetic transition as measured while warming 
after field cooling to 5 K. 

Our results and interpretation are in qualitative agree- 
ment with a recent report by M. Pregelj et al£ on the 
similar Cu3Bi(Se0 3 )20 2 Br compound [T c = 27.4 K). 
Small differences in the magnetization saturization value 
with H || c is likely linked to defects arising from differ- 
ent growth conditions. We normalize our susceptibility 
curves to moles of Cu while M. Pregelj et al. normalize to 
moles of formula unit (scaling factor of three difference) . 

The origin of the anomaly in l/x(T), first observed 



in powder samples by P. Millet et al.£ and subsequently 
observed by us along multiple axes of single crystals of 
Cu 3 Bi(Se03)20 2 Cl (Fig.© and Cu 3 Bi(Se0 3 ) 2 2 Br£ re- 
mains an unresolved issue. Below 150 K, the plot of 
H/M exhibits a slight curvature. We therefore conclude 
that the data is no longer accurately described by the 
Curie-Weiss law below 150 K (a dashed line fit below 
150 K remains in Fig. [5] to exemplify the discrepancy). 
Such deviations from Curie- Weiss law can be expected 
for quasi-two-dimensional spin systems and for frustrated 
topologies where non-trivial spin-correlation functions 
characterize the "classical spin liquid regime" or "corre- 
lated paramagnet" between the Curie- Weiss temperature 
(from the Curie- Weiss fit above 150 K) and the magnetic 
ordering temperature (24 K). 

We observed the anomaly for external fields of 1 T and 
0.01 T in both zero field cooling as well as field cooling 
measurements. 




Temperature (K) Temperature (K) 



FIG. 8. (Color online) Plot of H/M as a function of temper- 
ature measured with a 1 T magnetic field oriented parallel to 
the a axis (left panel), and also parallel to the c axis (right 
panel) . 



IV. DISCUSSION 

A. Group Theory and Observed Modes 

The number of optical modes in the 
Cu 3 Bi(Se0 3 )202Cl compound can be determined 
by group theory analysis. Using the SMODES 14 
program, we arrive at the following distribution of 
modes: 

T opUcal = + UB { 2 l R ] + UB^ 

+12A^ + 6B[f + 9B { 2 f + 12B ( 3 f + 9A lu 

where (R) and (IR) denote respectively Raman active 
and infrared active modes. The B 3u , B2„, and Bi„ modes 
are infrared active along the a, b, and c crystal axes re- 
spectively. The 9Ai„ modes are silent. The 12A g modes 
are Raman active radial breathing modes and require 
parallel incoming and outgoing polarizations according 
to the selection rules of the Pmmn spacegroup. The B 3g , 
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B2 ff , and Bi g modes are also Raman active, but require 
crossed polarization of the light (6c, ac, and ab respec- 
tively) . 

At 300 K we observe all 11 of the 11 predicted B^ u 
modes along the a axis, all 14 of the 14 predicted B2« 
modes along the b axis, and 11 of the 14 predicted Bi„ 
modes along the c axis. The discrepancy between ob- 
served and detected modes along the c axis is likely an 
experimental shortcoming arising from the low signal to 
noise when measuring along the thin dimension of the 
crystal. At 115 K we detect 8 additional modes along 
the a axis, 6 additional modes along the b axis, and 2 ad- 
ditional modes along the c axis. To study the observed 
modes, the reader is referred to Table U where the modes 
at 7 K along all three crystal axes are identified as well 
as their respective Lorentz oscillator parameters. An as- 
terisk after the TO frequency indicates a new mode that 
arises at 115 K. 

In materials with an inversion center, the k = op- 
tical modes as measured through Raman and infrared 
reflectance are mutually exclusive. When the inversion 
center in a structure is removed, the local centers of sym- 
metry about which all the Raman modes have zero net 
dipole moment will be removed, and all Raman modes 
will therefore become infrared active^ To further in- 
vestigate the 16 new modes observed at 115 K, prelimi- 
nary 300 K Raman measurements were recorded^ in the 
ab plane of crystal 1. Raman modes observed at 172.9, 
323.5, and 484.4 cm -1 closely correspond to three of the 
new infrared modes observed below 115 K. If we take into 
account the fact that Raman modes can slightly increase 
their resonance frequencies upon cooling to low tempera- 
tures, then 300 K Raman modes at 538.3 and 583.1 cm -1 
draw close redolence to two additional new modes ob- 
served in the infrared below 115 K. We therefore turn 
to a closer examination of a possible centrosymmetric 
to non-centrosymmetric transition in Cu3Bi(Se03)202Cl 
near 115 K. 



B. Powder X-Ray diffraction 

An outstanding problem in the interpretation of our 
results is the fact that the powder x-ray diffraction mea- 
surements did not find any evidence for a phase transition 
accompanying the dramatic appearance of 16 infrared- 
active modes below 115 K. In their 2001 report, P. Mil- 
let et al. A sought to detect a potential structural phase 
transition by optical birefringence, and they placed an 
upper limit of 1° on the possible rotation of optical axis, 
thus indicating no distortion to monoclinic symmetry at 
that level. The present powder diffraction measurements 
would be sensitive to a monoclinic distortion on the order 
of 0.001°, and none is observed. Their analysis of linear 
birefringence suggest that the nature of the suspected 
transition is second-order, which is consistent with what 
we have found in our analysis of shifts in spectral weight 
from existing modes to new modes (inset lower panel 



Figj3]). In what follows, we consider lower-symmetry non- 
centrosymmetric orthorhombic structures that might be 
a potential host for Cu3Bi(Se03)202Cl. In addition, we 
scrutinize the possibility of a transition to an incommen- 
surate lattice, and we discuss the implications of the re- 
cent neutron diffraction report on Cu3Bi(Se03)202Br. 5 

The centrosymmetric Pmmn space group has extinc- 
tion class P--n, i.e., it obeys the condition that h + k 
must be even for (hkO) reflections. A continuous transi- 
tion to a non-centrosymmetric orthorhombic subgroup of 
the same lattice dimensions could lead to space groups 
P2i2i2, Pmm2, Pm2\n : or P2\mn. The two former 
choices would allow additional powder x-ray diffraction 
peaks. We have carefully searched for them throughout 
the temperature range below 115 K without success. The 
two latter space groups have the same extinction class as 
Pmmn, and so they represent a mechanism for breaking 
of inversion symmetry without producing a qualitative 
change in the powder diffraction pattern. Such a distor- 
tion might not be easy to recognize from powder diffrac- 
tion data because the atoms would presumably move a 
small distance from their undistorted locations, and so 
the Pmmn model would still give a reasonably accurate 
description of the data in the acentric phase. One pos- 
sibility is to look for unusual behavior in thermal dis- 
placement parameters, similar to the method used to de- 
code complicated distortions in magnetically frustrated 
spinels. 17 We have not been able to detect such an effect 
from the data at hand. 

Neutron powder diffraction measurements on 
Cu3Bi(Se03)202Br at 60 K were refined in Pmmn^ 
but that has only limited bearing on the present issue. 
First, it is not known that data from the bromide 
material suggests the same loss of inversion symmetry 
that we have reported here for the chloride material. 
Second, as noted above, the difference between centric 
and acentric structures could be very difficult to see 
in neutron or x-ray powder diffraction. Indeed, it is 
worth pointing out that the magnetic transition studied 
in the bromide analog may well have occurred from a 
symmetry lower than Pmmn£&- 

Another possibility is that an incommensurate transi- 
tion occurs in Cu3Bi(Se03)202Cl below 115 K. Incom- 
mensurate phases acquire satellite x-ray diffraction lines 
around the Bragg peaks of the symmetric phased We 
did not observe satellite peaks, but it is possible they 
exist and were too weak to be detected in the present 
powder x-ray diffraction measurements. We plan to per- 
form more sensitive measurements (e.g., single crystal or 
neutron powder diffraction) to resolve this issue. 



C. Phonon Repulsion 

The typical dispersion of phonon resonance frequencies 
(i.e., softening with increasing temperature) as dictated 
by the anharmonic term in the lattice potential is not ob- 
served for a number of the new modes and existing modes 
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TABLE I. Oscillator parameters for the infrared observed modes of Cu3Bi(Se03)202Cl (at 7 K) along all three crystal axes. 
The new modes arising below 115 K are indicated with an asterisk next to their corresponding TO frequencies. 
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in close proximity to the new modes appearing below 
115 K. According to our presumption that Raman modes 
become infrared active with a loss of inversion symme- 
try, we examined the physics behind phonons of similar 
strengths resonating at contiguous frequencies. Figure [9] 
depicts a new mode appearing around ^328 cm -1 along 
the a axis that is in close proximity to an existing mode at 
^323 cm -1 . As further emphasized in the inset of Fig.Hl 
the two modes strongly repel one another as temperature 
decreases, which seems uncharacteristic of bosonic excita- 
tions. The repulsion is due to phonon mixing and can be 
understood by an analogy to the classical system of two 
coupled ideal harmonic oscillators. The classical problem 
amounts to solving for the eigenvalues of a matrix with 
off diagonal terms arising from the coupling between the 
oscillators. If both oscillators are given the same initial 
frequency (analogous to both phonons resonating at the 
same energy) and the coupling is turned on, then the os- 
cillators will repel one another. The repulsion increases 
with increasing coupling. Applying the same concept to 
the case of phonon mixing, one can see that the cou- 
pling between the phonons is increased with decreasing 
temperature, an effect which is well understood via the 
thermal contraction of the lattice. It is also worthwhile 
noting that the new phonon mode depicted in Fig. [3] at 
^276 cm -1 causes a repulsion of the resonant frequency 
of the existing mode at ^256 cm -1 below 115 K. The two 
modes appear to be initially non-degenerate because of 
their separation in energy (~20 cm -1 ); however, as seen 
in the inset of Fig. [3] lower panel, there is certainly an 
interaction between the two modes. The coupling is ver- 



ified further by the shifts in oscillator strengths between 
the two modes (lower panel Fig. [3] inset). 
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FIG. 9. (Color online) The real part of the optical conductiv- 
ity along the a axis in the frequency range 300-340 cm -1 . At 
115 K a new mode appears around ~328 cm -1 that is in close 
proximity to an existing mode at ~323 cm - . The repulsion 
of the two modes as temperature is decreased is emphasized 
in the inset where the Drude-Lorentz resonance frequencies 
are plotted as a function of temperature. 



D. Magnetic Excitations 

Magnetic excitations stimulated by infrared light can 
be grouped into one of two categories: (a) traditional 
magnons that are excited by the a.c. magnetic field of 
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the light, and (b) novel electromagnons that are ex- 
cited by the a.c. electric field component of the light. 
Single-magnon absorptions, which at infrared frequencies 
are commonly observed in antiferromagnets and thusly 
tabbed antiferromagnetic resonance modes (AFMR) , are 
magnetic dipole transitions that occur when the oscil- 
lating frequency of the light corresponds to the k = 
frequency of a spin wave. Electromagnons were first pro- 
posed as strongly renormalized spin waves with dipolar 
momentum^ Since their discovery, electromagnons have 
been extensively studied in rare earth manganite com- 
pounds, namely in TbMnC>3, where a large body of recent 
literature (Raman^, neutron^, and infrared—) has tied 
the excitation to the lattice itself. The work has resulted 
in a generalized hybrid magnon-phonon mode picture of 
electromagnons. Challenges arise in differentiating the 
two aforementioned excitations because they resonate in 
the same general frequency intervals. A common solu- 
tion is to measure the different faces of a crystal while 
rotating the polarization of the incoming light. 

The geometry of the Cu 3 Bi(Se03)202Cl crystal only 
allowed for transmission measurements with the k vec- 
tor of the incoming light aligned perpendicular to the 
ab plane. However, a thorough polarization study 
within the ab plane was carried out and it has lead 
to the observation of an isotropic magnetic excitation 
at 33.1 cm -1 . These results contradict the previously 
reported anisotropic nature of magnons and electro- 
magnons. Moreover, when external magnetic fields are 
applied in the H || c and H _!_ c geometries, the 33.1 cm -1 
resonance shifts to higher and to lower frequencies respec- 
tively. In what follows we will examine the nature of the 
excitation as well as propose a reason for its isotropic 
behavior. External field dependent spectra with H \\ c 
and H _L c will be discussed separately. (The following 
section will focus on the 33.1 cm -1 excitation because it 
is the only mode observed in zero field.) 



1. Nature and Isotropy 

To elucidate the nature of the magnetic excitation ob- 
served in Cu 3 Bi(Se03)20 2 Cl at 33.1 cm -1 , we compare 
the strength of the observed excitation to the exten- 
sive literature on magnons and electromagnons in the 
infrared. Our magnetic excitation creates a peak in 
a(ui) that is ~ 30 cm -1 above the baseline. Taking 
the average d.c. index of refraction to be 3 in the ab 
plane, we determine that our magnetic excitation corre- 
sponds to an optical conductivity of about 0.23 f2~ 1 cm _1 
(<ti = £^na), which is roughly the same strength as 
other reported single-magnon excitations*^ — Electro- 
magnons, which have been extensively studied in a num- 
ber of rare-earth manganites, are observed to possess op- 
tical conductivities of at least a factor of 10 larger.^S 
Although this method of comparison does not yield ob- 
jective certainty, we can further support its conclusion 
by employing optical sum rule analysis on our measured 



reflectance data. Since an electromagnon contributes to 
the dielectric constant, it must gain spectral weight from 
a dipole active excitation, the main candidates being do- 
main relaxations, phonons, or electronic transitions. The 
reflectivity spectra measured at 7 K and 30 K do not 
show any change associated with the magnetic excita- 
tion. Therefore, the magnetic excitation does not gain 
spectral weight from the low frequency infrared-active 
phonon modes. 



E||a 




E||a 



FIG. 10. (Color online) A polar plot of the oscillator strengths 
associated with the 33.1 cm -1 mode at T for all four polar- 
izations measured in the ab plane. 

The isotropic character of the magnetic excitation ob- 
served at 33.1 cm -1 is depicted in Fig. [TU] where the os- 
cillator strength of the excitation is plotted versus the 
polarization angle of the light in the ab plane. (It is 
worthwhile noting that infrared reflectance spectra mea- 
sured on crystal 2 revealed the same strong anisotropy 
as depicted in Figfl] thus excluding twinning of the 
surface.) A recent report on TbMn03 by Pimenov et 
aZ.— details the observation of a magnon and an electro- 
magnon, active along perpendicular directions, resonat- 
ing at the same frequency. This occurrence gives the 
illusion of an isotropic magnetic excitation and is wor- 
thy of consideration in Cu3Bi(Se03)202Cl. Strongly op- 
posing this argument are the nearly equivalent oscillator 
strengths measured in any two orthogonal polarizations, 
as seen in Fig. [TU] Electromagnons typically have oscil- 
lator strengths that are at least one order of magnitude 
greater than traditional magnons, as noted in the com- 
parison of oscillator strengths above. 

A second more plausible explanation is the occur- 
rence of weak magnons at two orthogonal polariza- 
tions. N. Kida et a/— observed a similar phenomena 
in DyMnOs, namely, weak excitations arising with the 
a.c. magnetic field oriented along both the a and c crystal 
axes. They supported their interpretation of two orthog- 
onal magnons by inelastic neutron scattering experiments 
on a similar rare-earth maganite where it was reported 
that two magnon dispersions curves from orthogonal axes 
crossed k = at the same energy. In Cu3Bi(Se03)202Cl, 
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we suspect that magnon dispersion curves from the [100] 
and [010] directions cross k = at 33.1 cm -1 (~4 meV); 
however, inelastic neutron scattering measurements are 
needed to support our hypothesis. 



2. H |[ c field dependence 

To analyze the ff || c spectra (i.e., H parallel to the 
easy axis), we will utilize our d.c. susceptibility measure- 
ments as well as generalize the magnetic structure de- 
termined for Cu 3 Bi(Se03)20 2 Br to Gu3Bi(Se03) 2 2 Cl. 
At zero field, six distinct magnetic sublattices can be 
identified, which would presumably lead to six distinct 
magnon branches; however, due to the canted nature of 
the copper ions occupying the 4(c) sites, the exact num- 
ber of branches could differ from six. At zero field, we 
only resolved one magnetic excitation (33.1 cm -1 ). No 
clear signature of the metamagnetic transition can be 
identified when tracking this excitation at fields between 

and 1 T. In a 1 T field, Cu3Bi(Se0 3 )20 2 Cl has al- 
ready undergone a metamagnetic transition where mag- 
netic moments on every second layer flip, resulting in 
ferromagnetic interlayer and canted ferrimagnetic behav- 
ior overall^ (It should be noted that the c axis remains 
the easy axis after the transition). Ironically, although 
the metamagnetic transition effectively reduces the num- 
ber of magnetic sublattices from six to three, we observe 
an additional magnetic excitation appearing at fields of 

1 T and greater. At this point, two logical questions 
arise: First, why would the 33.1 cm -1 excitation present 
in the low field antiferromagnetically-ordered state per- 
sist smoothly through the transition state and into the 
high-field fcrro or ferrimagnetically ordered state? Sec- 
ond, why are more magnetic excitations present in the 
high field phase despite a reduction in the number of 
magnetic sublattices? 

As to the first question, magnons resulting from ferro- 
magnetic resonance have been observed at infrared fre- 
quencies; however, they typically resonate at much lower 
frequencies because they must extrapolate linearly to 
zero frequency at zero applied magnetic field. It is there- 
fore very unlikely to observe a ferromagnetic resonance 
at 33.1 cm -1 in a 1 T external field. But there is an ex- 
ceptional case that was first observed by Jacobs in FeCl2 
apropos metamagnetic transitions. 51 When the metam- 
agnetic transition occurs in FeCl2 , which signals a transi- 
tion from a two-sublattice antiferromagnet to a ferromag- 
net, the AFMR line disappears in favor of a ferromagnetic 
resonance line around the same frequency. The high fre- 
quency of the ferromagnetic resonance line in FeCl2 is 
explained by large anisotropy fields in the material, and 
theoretical calculations support its existence^ We sus- 
pect that a similar situation occurs in Cu3Bi(Se03)202Cl 
that explains the smooth movement of the excitation at 
33.1 cm from below to above the metamagnetic tran- 
sition. Likewise, we suspect that the magnetic excita- 
tion at 10.5 cm -1 (1 T) also originates from an antifer- 



romagnetic resonance line, which would possess a zero 
field resonance frequency of 9.5 cm" 1 (slightly below our 
measurable range). 

Apropos the second question, that is to say, the pres- 
ence of more excitations above the metamagnetic transi- 
tion versus below it contradicts the reduction in magnetic 
sublattices from six to three. We suspect that more ex- 
citations do exist beneath the metamagnetic transition 
that are either below our measurable frequency range or 
that are too weak for us to resolve (resolution 0.3 cm -1 ). 
We have thoroughly inspected the low frequency spectra 
as well as the 33.1 cm -1 excitation between and 1 T, 
and we observe no clear signature of new resonances or 
the splitting of the existing 33.1 cm -1 excitation. It is 
possible that the 10.5 cm -1 excitation (at 1 T) splits be- 
neath the metamagnetic transition; however, below 1 T 
experimental limitations prevented us from tracking the 
excitation. Future studies involving a high-resolution, 
low-frequency source are needed to investigate further. 

In the interval 1-10 T we can estimate the effective 
g factor from the slope of the observed resonance lines 
at 10.5 and 33.1 cm -1 using the formula hw = g[isH. 33 
The slopes of the resonance lines at 10.5 and 33.1 cm -1 
correspond respectively to g factors of 2.16 and 0.24. 



3. H _L c field dependence 

Again, we use our d.c. susceptibility measurements and 
the neutron scattering results on Cu 3 Bi(Se03) 2 2 Br to 
conclude that antifcrromagnetic interactions exist within 
the ab plane at all measurable fields in the H _!_ c ge- 
ometry. The theory of antiferromagnetic resonances for 
uniaxial or cubic antiferromagnetic crystals was worked 
out by Keffer and Kittelj^ who showed that when the 
external field is oriented perpendicular to the easy axis, 
the two k = magnon branches correspond to frequen- 
cies w/7 = ±[2HeHa + Hq] 1 / 2 . In the previous equa- 
tion, 7 = ge/2mc where g is the spectroscopic split- 
ting factor, and Hq, Ha, and He represent the static, 
anisotropy, and exchange fields respectively. Generaliz- 
ing this theory, which was worked out for two sublattices, 
to a situation with more than two sublattices, we can see 
that the dispersion of the excitation we observe (inset of 
Fig. [5)3) is in qualitative agreement (i.e., both are non- 
linear) with the lower frequency branch predicted by the 
AFMR theory. (Further experiments are needed to de- 
termine Ha , He , and g and subsequently fit the formula 
to the data to either validate or undermine our general- 
ization.) When there is no external field, the two absorp- 
tions are predicted to be degenerate. However, in antifer- 
romagnets with strong magnetic anisotropy, which is the 
case for Cu 3 Bi(Se03)202Cl, the degeneracy of magnon 
branches in zero field is lifted (e.g., MnO^ NiO^ and 
NiF2^). Following our previous reasoning that the res- 
onance observed in Cu3Bi(SeO-3)202Cl, which moves to 
lower frequency with increasing field, is the lower fre- 
quency branch of the k = AFMR, we conclude that the 
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higher frequency branch becomes masked, thus unobserv- 
able, behind the strong phonon absorptions starting ^40 
cm -1 . AFMR theory also predicts that when the lower 
frequency branch reaches zero the external field is able to 
rotate the spins, at the expense of the anisotropy field, 
away from the easy axis (c axis) and into the ab planed 
Extrapolating the lower frequency branch towards zero, 
while taking into account the uncertainty in the data 
points, we estimate a zero frequency crossing somewhere 
between 18.5 and 20.1 T. Our predicted field range at 
which the spins rotate away from the easy axis is slightly 
higher, but remains in rough agreement with the values 
predicted on the Cu3Bi(Se03)202Br analogue^ But un- 
like Cu3Bi(Se03)202Br, we can not identify intermedi- 
ate and hard axes because the movement of our observed 
mode does not seem to depend on the orientation of the 
external field within the ab plane. 

V. CONCLUSIONS 

The novel geometrically-frustrated layered compound 
Cu3Bi(Se03)202Cl has been characterized using infrared 
spectroscopy, powder x-ray diffraction, and d.c. mag- 
netic susceptibility measurements. Far-infrared re- 
flectance measurements have revealed 16 new infrared- 
active phonon modes below 115 K. The plethora of 
new modes observed strongly suggest a rearrangement of 
atomic positions within the unit cell; however, our sub- 
sequent powder x-ray diffraction measurements are com- 
pletely consistent with the same 300 K structure (Pmmn) 
existing at 85 K. Preliminary Raman spectra taken at 
300 K on crystal 1 have revealed five phonon modes at 
frequencies close to five of the new modes observed in 
the infrared below 115 KJ£ The results suggest a loss 
of inversion symmetry below 115 K. Upon further inves- 
tigation we have identified two non-centrosymmetric or- 
thorhombic space groups {Pm2\n and P1\mn) that have 
the same allowed Bragg reflection peaks as the 300 K 
Pmmn structure. Therefore we suspect that a subtle 
second-order transition from Pmmn to either Pmflin or 
P2\mn occurs near 115 K that is below the resolution 



of our powder x-ray diffraction experiment. We plan to 
perform more sensitive measurements (e.g., single crystal 
or neutron powder diffraction) to resolve the uncertainty 
of this issue. 

In addition, an isotropic magnetic excitation is ob- 
served at 33.1 cm -1 at 5 K. We have tentatively assigned 
the magnetic excitation to a magnon based on analysis of 
previously reported oscillator strengths of magnons and 
electromagnons. The isotropic behavior of the excitation 
within the ab plane is potentially due to the Brillouin 
zone center crossing of two magnon dispersion curves 
along orthogonal directions, but inelastic neutron scat- 
tering measurements are needed to investigate further 
the excitation's seemingly isotropic existence. 

The resonance frequency of the 33.1 cm -1 excita- 
tion strongly depends on the orientation of the static 
magnetic field. An anisotropic response to the orienta- 
tion of a static magnetic field is also seen in d.c. sus- 
ceptibility measurements on Cu3Bi(Se03)202Cl, as well 
as neutron diffraction measurements on the similar 
Cu3Bi(Se03)202Br compound. 

When the external magnetic field is applied parallel 
to the c axis (H || c), the resonant frequency of the 
33.1 cm -1 excitation increases linearly with increasing 
field. For fields of 1 T and greater applied along the c 
axis an additional linearly-increasing magnetic excitation 
is observed (10.5 cm -1 at 1 T). 

When the external magnetic field is applied perpen- 
dicular to the c axis (H _L c), the resonant frequency 
of the 33.1 cm -1 excitation decreases quadratically with 
increasing field. The results are in agreement with the 
behavior of an antiferromagnetic resonance line in the 
presence of strong magnetic anisotropy. 
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